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Abstract. This paper reviews the most recent developments related 
to the use of lensing clusters of galaxies as Gravitational Telescopes in 
deep Universe studies. We summarize the state of the art and the most 
CO . recent results aiming at studying the physical properties of distant galax- 

ies beyond the limits of conventional spectroscopy. The application of 
■ photometric redshift techniques in the context of gravitational lensing is 

Q-f emphasized for the study of both lensing structures and the background 

q ■ population of lensed galaxies. A presently ongoing search for the first 

\ building blocks of galaxies behind lensing clusters is presented and dis- 

cussed. 

03 
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. 1. Introduction 

The availability of statistically significant samples of galaxies, from z ~ to 
the largest look-back times, covering a large region in the parameter space, 
is mandatory for constraining the cosmological scenarios of galaxy formation 
and evolution. Large spectroscopic samples of galaxies at all redshifts have 
become available during the last ten years, thanks to extended surveys in the 
different rest-frame wavelength domains (e.g. Lilly et al. 1995; Ellis et al. 1996; 
Steidel et al. 1996, 1999; Adelberger & Steidel 2000; Chapman et al. 2000; 
Pettini et al. 2001; Erb et al. 2003). The evolution of the overall properties of 
galaxies as a function of redshift is particularly important for galaxies at z > 1, a 
redshift domain where galaxies are expected to be strongly affected by merging 
or assembly processes. 

Clusters of galaxies acting as Gravitational Telescopes (hereafter GTs) con- 
stitute a particular and powerful tool in this context of high redshift surveys, as 
it was pointed out by different authors in the early nineties (see Fort &: Mellier 
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Figure 1. Spectra of three emblematic lensed LBGs identified in the 
nineties, showing the most prominent absorption lines. From top to 
bottom: #384 in A22f8 (z= 2.515, Ebbels et al. 1996), S2-AC114 
(z= 1.867, Smail et al. 1995; Campusano et al. 2001), and cB58- 
MS1512+36 (z=2.7, Yee et al. 1996; Pettini et al. 2000). 



1994 for a review). The main advantage of such "telescopes" is that they take 
benefit from the large magnification factor in the core of lensing clusters, close 
to the critical lines, which typically ranges between 1 and 3 magnitudes. Thus, 
GTs can be used to access the most distant and faintest population of galaxies, 
extending the current surveys towards intrinsically fainter samples of galaxies, 
allowing to probe the physical properties of the faint end of the luminosity func- 
tion at moderate and high redshifts. Indeed, GTs have been succesfully used 
during the last ten years to perform detailed studies of distant galaxies at differ- 
ent wavelengths, from UV to submillimeter. A particular and recent application 
of GTs is the identification and study of extremely distant galaxies, including 
the search for the first galaxies formed in the early Universe, up to redshifts of 
the order of z ~ 10 — 20. 

The plan of the paper is as follows. Section [2D presents a summary of the 
latest developments concerning the use of cluster lenses as Gravitational Tele- 
scopes. We summarize the state of the art and the recent results aiming at 
studying the physical properties of distant galaxies beyond the limits of conven- 
tional spectroscopy In Section |23 we briefly describe the photometric redshift 
approach, as a tool in deep photometric studies, and in particular in the context 
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Figure 2. Velocity field in AC114-S2 (z= 1.867, from Lemoine- 
Busserolle et al. 2003). Lens corrected velocity profile (in km/s) versus 
distance measured on the 2D spectra from the central velocities of Ha 
(open dots) and [OIII] A5007 (full dots) emission lines. A velocity 
gradient of ±240 km/s is measured, which yields a dynamical mass of 
~ 1.3 x 10 10 Mq within the inner 1 kpc radius. 



of gravitational lensing. We emphasize the use of such a tool for the study of 
both lensing structures and the background population of lensed galaxies. We 
discuss in Section ^\ on the use of Gravitational Telescopes to identify the first 
building blocks of galaxies, also known as Population III sources. A presently 
ongoing search for PopIII candidates behind lensing clusters is briefly presented. 
Conclusions are given in Section HT1 



2. Physics of galaxies beyond the limits of conventional spectroscopic 
samples. State of the art 

2.1. Lensing Clusters as Gravitational Telescopes 

From the point of view of galaxy evolution studies, GTs can be considered as an 
alternative way to blank fields to investigate the properties of distant galaxies. 
GTs allow to construct and to study independent samples of high-z galaxies, less 
biased in luminosity than standard magnitude-limited field samples, thanks to 
the magnification factor. For a given limiting magnitude or flux, lensed samples 
will complement the currently available field surveys towards the faint end of 
the luminosity function, and towards higher limits in redshift. It is worth noting 
that only lensing clusters with fairly well constrained mass distributions can 
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be used as efficient GTs for galaxy evolution studies. In order to optimize the 
surveys, the expected 2D surface density distribution of arclets in clusters can 
be estimated for lenses with well known mass-distributions. Conversely, the 
lens-corrected distribution of arclets can be easily retrieved from the cluster 
cross section (e.g. Kneib et al. 1994, Bezecourt et al. 1998). An interesting 
property of GTs is that they conserve the surface brightness and the spectral 
energy distribution (SED) of lensed galaxies. Besides the spatial magnification, 
sources (or regions of sources) smaller than the seeing will effectively gain in 
surface brightness. These two effects contribute to increase the interest of GTs 
in spectroscopic studies. 

The use of lensing clusters as genuine GTs is a relatively new research 
domain. Historically, it started with the spectroscopic confirmation of the first 
highly amplified source at z = 0.724 in A370 (Soucail et al. 1988), and the 
identification of the spectacular blue arc in C12244-02 as a high-z lensed galaxy at 
z ~ 2 (Mellier et al. 1991). The power of GTs to study the spectrophotometrical 
and morphological properties of high-z galaxies was emphasized by different 
authors in the nineties (e.g. Smail et al. 1993; Smail et al. 1996; Kneib et al. 
1996). 

One of the first detailed studies on the physical properties of high-z lensed 
galaxies was carried out by Ebbels et al. (1996), who studied the optical spec- 
trum (obtained with the 4m telescope WHT) and optical-infrared colours of a 
star-forming galaxy at z=2.515 behind A2218. This galaxy exhibits the same 
properties as typical Lyman Break Galaxies (LBGs, Steidel et al. 1996) in terms 
of colors, luminosity and star formation rate (see also Leitherer 1999). Another 
emblematic source is the luminous fold arc cB58, a z=2.7 star-forming galaxy 
behind the EMSS cluster MS1512+36, which is probably the best studied LBG 
nowadays (Yee et al. 1996; Seitz et al. 1998; Pettini et al. 2000; Teplitz al. 
2000). Four galaxies with 3.3 < z < 4 were found in the rich cluster C10939+47 
by Trager et al. (1997), three of them at z ~ 4. Because of the weakness of 
the common stellar wind features in the spectra of these objects, the authors 
concluded that they have subsolar metallicities, their properties beeing quite 
similar to the already known ones in LBGs at z ~ 3. Because these objects, 
revealed by the GT, are close together both in redshift and space, they were 
identified by the authors as low-metallicity Population II protospheroid clumps, 
in other words, the progenitors of today's spheroids. Another close look towards 
a star-forming region at z ~ 4 was provided by the two lensed galaxies found in 
the cluster Abell 2390, at z = 4.05 (Frye & Broadhurst 1998; Pello et al. 1999). 
These galaxies are relatively bright compared to LBGs at z > 3, or to the z ~ 4 
sources in C10939+47. They have subsolar metallicities according to the overall 
fit of their SEDs, and they are clumpy and elongated, as revealed by the lens 
inversion procedure. They are located ~ lOO/igQ 1 kpc apart from each other on 
the source plane, with similar orientation. Again, these results point towards 
a merging process of protospheroidal clumps, revealed by the GT. One of the 
highest redshift objects known in the nineties was also discovered thanks to the 
GT approach. This object was found in the cluster CL1358+62, with z = 4.92 
(Franx et al. 1997). Its spectral energy distribution is quite blue, as expected 
for a young galaxy with the observed restframe UV flux, but Soifer et al. (1998) 
results pointed towards an important absorption by dust in this object, the first 
direct evidence for such effect at z ~ 5. 
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Figure 3. Oxygen abundance as a function of Mb, showing the loca- 
tion of the two z ~ 1.9 lensed galaxies S2 and A2 in AC114 (from 
Lemoine-Busserolle et al. 2003), with (filled squares) and without 
(empty squares) extinction correction. The metallicity-luminosity re- 
lations for nearby "normal" (solid line; Kobulnicky & Zaritsky 1999) 
and starburst (dashed line; Mouhcine & Contini 2002; Melbourne & 
Salzer 2002) galaxies are also shown, together with different sam- 
ples of intermediate-redshift galaxies. Blue squares are high-redshift 
(1.4 < z < 3.4) galaxies with O/H derived from rest-frame UV spectra 
by Mehlert et al. (2002). LBGs at z ~ 3 are shown as a blue box 
encompassing the range of O/H and Mg derived for these objects by 
Pettini et al. (2001). 
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GTs have been also used to address the morphological properties of high-z 
lensed galaxies, in particular using HST images after its refurbishment. Grav- 
itational arcs on these images are detected with diffraction-limited resolution. 
Therefore, morphological criteria have been introduced to identify lensed sources 
(Kneib et al. 1996). Smail et al. (1996) used this approach to discuss the prop- 
erties of a large number of arcs in well known clusters, and they concluded to 
a strong evolution in the size of star-forming galaxies between z ~ 1 and the 
local universe, distant sources beeing 1.5 to 2 times smaller than local ones. A 
detailed multi- wavelength study of a ring- like starburst galaxy at z ~ 1, hosting 
a central active nucleus, was presented by Soucail et al. (1999) as a result of 
their survey in the field of A370. The comparison between this galaxy and other 
similar galaxies in the local universe was possible thanks to lens inversion. Also, 
as mentioned above, several studies of the close environment of z ~ 4 lensed 
sources were conducted with GTs. 

When lensed galaxies are selected for subsequent spectroscopic studies, they 
are usually chosen close to the high-z critical lines in order to obtain the largest 
magnification (typically ~ 2 magnitudes). Number of high-z lensed galaxies ob- 
served during the nineties were discovered serendipitously, but efficient selection 
procedures can be defined, based on photometric or lensing criteria. Therefore, 
number of recent surveys have started using GTs to constitute large samples of 
intrisically faint objects at intermediate and high redshifts, using the capabil- 
ities of present day multi-object spectroscopy facilities in the optical domain. 
Among them, Campusano et al. (2001) presented a spectroscopic survey of 
10 faint lensed galaxies in the core of the cluster AC114 (z c / US te r =0.312), with 
0.7<z<3.5. Of particular interest are five objects (A, B, C, S and E) with 
spectroscopic redshift 2.0<z<3.5, which are between 0.5 and 1.5 magnitudes 
fainter than the limiting magnitude in the Steidel et al. (1999) sample of LBGs 
at similar redshifts. In a recent paper, Mehlert et al. (2002) discussed on the 
evidence of chemical evolution in the spectra of high-z galaxies, by comparing a 
local sample of galaxies with a VLT/FORS spectroscopic sample at 1.4 < z < 3.4, 
using the equivalent width of CIV as a metallicity indicator. Part of their sam- 
ple is constituted by high-z lensed galaxies observed in the field of 1E0657-558 
(Mehlert et al. 2001), with very large magnifications (~ 3 mags.) allowing 
to compare the abundances measured in faint galaxies in the local and high- 
z universe. These authors conclude that the major trend in their sample is a 
metallicity enhancement towards low-z galaxies, with no overall dependence of 
metallicity with galaxy luminosity. These results have been derived from a small 
number of intrinsically faint galaxies, and they have to be confirmed by larger 
samples of lensed galaxies at high-z. 

2.2. A few words about sample selection criteria 

During the last ten years, samples of high-z lensed galaxies have been pre- 
selected in the core of lensing clusters using different techniques. The fraction 
of high-z galaxies in magnitude-limited samples is relatively small, even in deep 
photometric surveys, but the same photometric selection techniques usually ap- 
plied to blank field studies can be adopted to select high-z galaxies behind lensing 
clusters. Broadband color criteria have been successfully used to identify galax- 
ies at redshifts higher than z ~ 2 (e.g. Steidel et al 1995; Erb et al. 2003), 
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Figure 4. Same caption as in Figured Oxygen abundance as a func- 
tion of line-of-sight velocity dispersion (a), from Lemoine-Busserolle et 
al. (2003). The dashed line is a linear fit to nearby and intermediate-z 
samples. LBGs at z ~ 3 are shown as a blue box encompassing the 
range of O/H and a derived for these objects by Pettini et al. (2001). 



and more general photometric redshift techniques are likely to produce equally 
usefull results in this context (see Section [TJ). In the paper by Campusano et 
al. (2001), lensed galaxies were selected in areas close to the high-z critical 
lines predicted by the gravitational lens model of Natarajan et al (1998), using 
a combination of both lensing and photometric redshift criteria. All galaxies 
targetted by the spectroscopic follow-up were found to be background galaxies 
with redshifts values within the 0.7 < z < 3.5 interval. 

Another selection criterium was proposed by Kneib et. al. (1994) and Kneib 
et. al. (1996) based on the inversion of the lensing equations and the prediction 
of redshifts for individual arclets. This is a purely geometrical criterium, based 
on the averaged shape of lensed sources, and the accurate measurement of the 
shape and orientation of arclets on high-resolution HST images, although the 
first applications were obtained on ground-based data (Kneib et. al. 1994). This 
technique produces statistical results on the redshift distribution of background 
sources. The spectroscopic verification of this technique on statistical basis was 
obtained by Ebbels et al. (1998) in the field of A2218. Bezecourt et al. (1999) 
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used this lens inversion method to derive the comoving star formation rate in 
the range 0.5<z<2 from UV (WFPC2/F336W) images of the cluster lens A370. 

2.3. The near-IR window 

Observations in the near-IR window have become widespread during the last 
years. Indeed, without this wavelength domain, only star forming galaxies will 
be entering the high-z samples of galaxies defined in the visible bands. Deep 
observations at longer wavelengths are needed to probe stellar masses, thus 
constraining the cosmic star formation history. The near-IR window is barely 
affected by the presence of dust extinction or starbursts, and k-corrections are 
reduced compared to optical bands. This window is required to consistently 
follow the stellar population contributing to the flux beyond A >4000 A from 
z ~ all the way to redshift 4, and to map the star-forming and AGN activity 
up to the highest redshifts (z ~ 10). The evolution of the overall near-IR proper- 
ties of galaxies (Luminosity Functions, redshift distribution of near-IR samples, 
color distributions, ...) as a function of redshift constitutes a powerful test to 
discriminate between the different scenarios of galaxy formation, in particular 
for galaxies at 1 < z < 2, a redshift domain where GTs are particularly efficient. 
The Hubble Deep Field data have shown that the fraction of galaxies at z > 1 
reaches 50% typically at K~22.5-23. This is indeed a key redshift domain where 
most galaxies were affected by major merging/assembly processes. In the case of 
early type galaxies, the formation of massive objects seems to took place at red- 
shifts z > 2, i.e. the majority of their masses were already assembled and they 
only underwent minor mergers at lower redshifts. The main spectral signatures 
of such assembly process could be identified and studied through optical and 
near-IR data, the later beeing crucial to determine the photometric redshifts of 
galaxies within the redshift interval 1 < z < 2 (see Section • 

The first results on a recent photometric survey for Extremely Red Objects 
(EROs), undertaken with the GT on a sample of 10 clusters at z ~ 0.2, was 
recently published by Smith et al. (2002). 60 EROs were found with (R — K) > 
5.3 down to K = 20.6. After correction for lensing effects, number counts are 
found to flatten for magnitudes fainter than K ~ 19-20, a result which could be 
due to the transition from a population of EROs dominated by evolved galaxies 
at z ~ 1-2 (K < 19-20) to one dominated by dusty starbursts at z > 1 (K > 19- 
20). 

The recent development of near-IR spectrographs on 10m class telescopes 
has allowed the study of the rest-frame optical properties of high redshift galax- 
ies. Indeed, near-IR spectroscopy allows to access the most relevant emission 
lines (Ha, H/3, [Oil] AA3726,3728, [OIII] AA4959,5007, ...), in order to probe the 
physical properties of galaxies (SFR, reddening, metallicity, kinematics, virial 
mass, etc), all the way from the local universe to z ~ 4, using the same parameter 
space and indicators. The pioneering work by Pettini and collaborators (1998, 
2001) has shown that the rest-frame optical properties of the brightest Lyman 
break galaxies (LBGs) at z ~ 3 are relatively uniform. A new sample of 16 LBGs 
at 2 < z < 2.6 has been recently published by Erb et al. (2003) and, in this case, 
significant differences are found in the kinematics of galaxies at z ~ 2 compared 
to z ~ 3. In a recent paper, Lemoine-Busserolle et al. (2003) have presented 
the firsts results on their near-IR spectroscopic survey of high-redshift magnified 
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galaxies at z ~ 2 using GTs. Two 2 ~ 1.9 lensed star- forming galaxies have been 
studied in the field of the cluster AC 114 (S2 and A2), 1 to 2 magnitudes fainter 
than LBGs at z ~ 3. Thanks to the large magnification factors obtained for 
these objects, number of rest-frame optical emission lines have been measured 
from [OII]A3727 to Ha+[N II]A6584. Before this paper by Lemoine-Busserolle 
et al. (2003), there was only one LBG for which chemical abundances have been 
determined with some degree of confidence: the gravitationally lensed galaxy 
MS 1512-cB58 at z = 2.729 (Teplitz et al. 2000; Pettini et al. 2002a). Figure 01 
displays the metallicity-luminosity relationship and the metallicity versus line- 
of-sight velocity dispersion (cr, which corresponds to a crude mass-metallicity 
sequence) for these lensed galaxies compared to reference samples of galaxies at 
different redshifts. The results obtained up to now on this dramatically small 
sample suggest that high- 2: objects of different luminosities could have quite dif- 
ferent histories of star formation. The sample of LBGs observed in the near-IR 
aiming to explore the physical properties of galaxies is still poor. The contribu- 
tion of GTs to these detailed studies is already significant, and there is a good 
chance that it increases in the following years. 

2.4. Extending the wavelength domain 

In recent years, GT studies of lensed galaxies have successfully extended to other 
wavelength windows: 

• In the mid- infrared (MIR), ultra-deep observations through the lensing 
cluster A2390 were obtained with ISOC AM/ISO at 7 and 15 n m (Altieri 
et al. 1999). A large number of sources were identified as high-z lensed 
galaxies, and number counts computed to unprecedent depth. These re- 
sults ruled out non-evolutionary models and favour of a very strong evo- 
lution (see also Metcalfe, L., this Conference). 

• In the sub-millimeter domain, the pioneering work by Smail et al. (1997) 
with SCUBA/JCMT, on the lensing clusters A370 and C12244-02, showed 
that the majority of lensed sources detected at 450 and 850 \x m lie at red- 
shifts above z ~ 1. In more recent papers (Smail et al. 1998; Blain et al. 
1998, 1999; Ivison et al. 2000; Smail et al. 2002; see also Kneib, this Con- 
ference), GTs have been used to trace the dust obscured star- formation ac- 
tivity at high-z, thanks to enlarged samples of faint submillimeter-selected 
sources. It appears from these developments that the bulk of the 850/iin 
background radiation could be originated in distant ultraluminous galax- 
ies, and that a large fraction of z < 5 sources could be missing from optical 
surveys. 

2.5. Identification and study of extremely high-z sources 

Nowadays, one of the most exciting challenges for GTs is the identification and 
study of extremely high-z sources. Number of devoted surveys, using different 
techniques, are specifically tuned up to search for z>5 galaxies using GTs (see 
also Kneib, this Conference). 

The first results of these systematic searches for z > 5 galaxies in GTs have 
appeared during the last two years. A z= 5.576 double-imaged galaxy has been 
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detected in the field of A2218 (Ellis et al. 2001). This source corresponds to an 
extremely faint and compact object (I~ 30), highly magnified (about a factor 
of ~ 30), with a SED dominated by the Lya emission line. According to the 
authors, this source could be a ~ 10 6 M Q stellar system undergoing its first 
generation of stars. Another galaxy at z= 6.56 has been identified in the field 
of the cluster A370 by Hu et al. (2002). In this case, the magnification is a 
factor of ~ 4.5 for this source, which sets strong constraints on the reionization 
epoch. We present in Section \^ a particular ongoing project devoted to the 
identification of Population III candidates in lensing clusters. 

3. Photometric Redshifts and Gravitational Telescopes 

3.1. Photometric Redshifts, a unique tool in deep universe studies 

Today, photometric redshifts (hereafter z p h ot s) constitute a well established tool 
in deep universe studies. This technique was originally proposed by Baum 
(1962), to measure the redshifts of elliptical galaxies in distant clusters, and 
rediscovered later by several authors in the eighties (Couch et al. 1983, Koo 
1985), this time applied to relatively low-redshift samples of galaxies observed 
in the ~ 4000 to 8000 A domain. The interest for this technique has increased 
later in the nineties, with the development of extremely large and/or deep field 
surveys. 

Indeed, deep photometric galaxy samples have become available during the 
last decade, in particular the observations by HST of the HDF-N (Williams et al. 
1996) and HDF-S (Casertano et al. 2000), and the coordinated complementary 
observations from the ground at near-IR wavelengths (Dickinson et al. 2001; 
Labbe et al. 2003). At the same time, reliable photometric redshift techniques 
have been developed by different authors, allowing distance estimates of faint 
galaxies for which no spectroscopic redshifts can be obtained nowadays, even 
with the most powerful telescopes (e.g. Connolly et al. 1997; Wang et al. 1998; 
Giallongo et al. 1998; Fernandez-Soto et al. 1999; Arnouts et al. 1999; Bolzonella 
et al. 2000; Furusawa et al. 2000; Rodighiero et al. 2001; Rudnick et al. 2001; 
Le Borgne &: Rocca-Volmerange 2002; and the references therein). 

3.2. Methods and accuracy 

There are basically two different techniques for measuring z p h ts : the so-called 
empirical training set method, and the SED fitting method. The first approach 
was originally proposed by Connolly et al. (1995, 1997). It consists on deriving 
empirically the relationship between magnitudes and redshifts using a subsample 
of objects with measured spectroscopic redshifts, i.e. the training set. A slightly 
modified version of this method was proposed by Wang et al. (1998) to derive 
redshifts in the HDF-N by means of a linear function of colours. The main 
advantages of this method are the small dispersion attained, even with a small 
number of filters, and the fact that it does not make use of any assumption 
concerning galaxy spectra or galaxy evolution, thus bypassing the problem of 
a poor knowledge of high redshift spectra. The main disadvantage is that the 
empirical relation between magnitudes and redshifts has to be recomputed for 
each filter set and survey, on a well suited spectroscopic subsample, thus making 
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this approach not very flexible. In addition, the redshift range between 1.4 and 
2.2 had been hardly reached by spectroscopy up to now, because of the lack of 
strong spectral features accessible to optical spectrographs. And efficient multi- 
object near-IR spectrographs are still missing on 10m class telescopes. Thus, no 
reliable empirical relation can be obtained in this interval. 

In the SED fitting procedure, the photometric redshift of a given object 
simply corresponds to the best fit of its photometric SED by the set of template 
spectra, obtained through a standard minimization procedure. This method ob- 
viously relies on the fit of the overall shape of spectra and the detection of strong 
spectral features, such as the 4000 A break, the Lyman decrement or particularly 
strong emission lines. In the nineties, this method has been applied to the HDF 
data, using either observed or synthetic template spectra (e.g. Mobasher et al. 
1996; Lanzetta et al. 1996; Gwyn & Hartwick 1996; Sawicki et al. 1997; Gial- 
longo et al. 1998; Fernandez-Soto et al. 1999; Arnouts et al. 1999; Furusawa et 
al. 2000; Bolzonella et al. 2000; Rudnick et al. 2001). A combination of this 
method with the Bayesian marginalization introducing a prior probability was 
successfully proposed by Bemtez (2000) (see also Coe, this Conference). The 
advantage of the SED fitting procedure is that no (large) spectroscopic training 
set is needed, although a control set is highly recommended. 

In order to estimate the accuracy of z p hots , a crucial test is the compar- 
ison between the photometric and the spectroscopic redshifts obtained on a 
restricted subsample of relatively bright objects, even if there is no guarantee 
that we are dealing with the same type of objects in the photometric (faint) and 
the spectroscopic (bright) samples. A complete discussion on the expected ac- 
curacy of Zphots , based on simulations and available spectroscopy in the HDFs, 
was presented by Bolzonella et al. (2000). In this paper the authors used the 
public code hyperz, which adopts a standard SED fitting method, but most re- 
sults should be generalized to other z p h ot tools, hyperz is available on the web 
at http : / / webast . ast . obs-mip . f r/hyperz. Figure El displays a comparison 
between spectroscopic and photometric redshifts in the HDF fields. The rms 
uncertainties on z p hot are typically a z ~ 0.05(1 + z), and up to a z ~ 0.1(1 + z), 
depending on the photometric depth, the redshift of the sources and the spectral 
coverage of the filter set. Without near-IR data, errors on z p hot are relatively 
higher in the redshift range 1 < z < 2, due to the lack of strong signatures in the 
observed SED between 3500A and 9000A. The 4000A break span the near-IR 
domain at redshifts between 1.2 and 5. The gain using near-IR photometry is 
also sensible for the determination of rough spectral types, such as early versus 
late type galaxies, or star-galaxy-qso discrimination. A complete discussion can 
be found in Bolzonella et al. (2000). 

3.3. Photometric Redshifts and lensed sources 

Photometric Redshifts techniques in lensing clusters allow to determine the red- 
shift distribution and properties of high-z and/or faint samples of lensed galaxies, 
sometimes invisible otherwise. GTs with well constrained mass distributions en- 
able to recover precisely the lens-corrected redshift distribution of background 
sources. Thus, z p h ot s allow not only to select high-z lensed sources for subsequent 
spectroscopic studies (as seen in Section[T]), but also to obtain robust constraints 
on the statistical properties of the background population of galaxies. 
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Figure 5. HST image of the cluster Abell 2390 in I band (F812W), 
showing the identification of the double-imaged source H5 at z= 4.05 
(from Pello et al. 1999). a) Mean spectra of the a and b components 
of H5, each of them showing a strong Lya emission line at z = 4.05.. 
b) Photometric redshift likelihood-map of H5 showing the excellent 
agreement with the spectroscopic redshift, which is contained within 
the region at 68% confidence level. The shaded region on the lower 
part of the Age vs. redshift map is excluded because of age-limit con- 
siderations for the stellar population (stars cannot be older than the 
age of the universe, with Hq = 50 km s _1 Mpc -1 and qo = 0). c) 
Restored I\y image of H5 on the source plane at z = 4.05, as obtained 
from the lens-inversion procedure. 
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One of the main issues is the determination of deep number counts, which 
are strongly related to the evolution of the spectral content of galaxies, a pro- 
cess driven by the star formation history. The cumulative redshift distributions, 
specially in the near-IR bands, can be used as a direct test for the scenarios 
of galaxy formation (Kauffmann Sz Chariot 1998; Fontana et al. 1999; Cole et 
al. 2000). Photometric Redshifts techniques have been recently applied to the 
study of galaxy properties up to the faintest magnitudes on the HDFs, to deter- 
mine the star formation history at high redshift from the UV luminosity density 
(Lanzetta et al. 2002), the evolution of the luminosity functions (Bolzonella et 
al. 2002), and also to analyse the stellar population and the evolutionary prop- 
erties of distant galaxies (e.g. SubbaRao et al. 1996; Gwyn &: Hartwick 1996; 
Sawicki et al. 1997; Connolly et al. 1997; Pascarelle et al. 1998; Giallongo et al. 
1998; Fernandez-Soto et al. 1999; Poli et al. 2001), or to derive the evolution of 
the clustering properties (Arnouts et al. 1999; Magliocchetti & Maddox 1999, 
Arnouts et al. 2002). All these developments can be easily applied to lensing 
clusters for which deep multicolor photometric surveys are available. 

A Monte-Carlo approach has been recently adopted by different authors to 
compute statistical quantities from photometric redshifts, such as the evolution 
of luminosity functions in deep photometric surveys, both in blank and strong 
lensing fields (SubbaRao et al. 1996; Dye et al. 2001). An interesting method 
was proposed by Bolzonella et al. (2002) to account for the non-gaussianity of 
the redshift probability functions, and specially to include degenerate solutions 
in redshift. 

An interesting issue for z p hotS associated to GTs for the spectroscopic follow- 
up of faint amplified sources is the optimization of the survey, i.e. selecting the 
best spectral domain in the visible or near-IR bands in order to observe the 
strongest spectral features. An additional benefit of z p h ot s is that this technique 
efficiently contributes to the identification of objects with ambiguous spectral 
features, such as isolated emission lines. An example of this is given in Campu- 
sano et al. (2001). 

An important concern for the different z p hot methods is that the training 
set (or the control set in the case of SED fitting) is constituted by the brightest 
objects, for which a spectroscopic measure of the redshift could be obtained. 
There is no guarantee that we are dealing with the same type of objects in the 
photometric (faint) and the spectroscopic (bright) samples. Thanks to the mag- 
nification factor, cluster lenses allow the calibration of z p hotS beyond the spec- 
troscopic limits in standard surveys, i.e. enlarging the training/control samples 
towards the faintest limits in magnitude. 

Another interesting and still poorly explored domain is the combination be- 
tween photometric and lensing redshifts. Indeed, lensing inversion and z p hot tech- 
niques produce independent probability distributions for the redshift of ampli- 
fied sources. Therefore, the combination of both methods shall provide a ro- 
bust way to determine the redshift distribution of the faintest sources, at least 
at z<l.5 — 2, where the lensing inversion technique produces reliable results 
(Ebbels et al. 1998). In general, the z p hot determination is more accurate than 
the zi ens for individual objects, in particular at z> 1.5. Nevertheless, the combi- 
nation of both probability distributions is particularly useful when a degenerate 
solution appears in z p hot ■ 
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3.4. Photometric Redshifts and lensing clusters 

Cluster lenses become useful GTs when their mass distribution is highly con- 
strained by multiple images, revealed by high resolution HST images or multi- 
color photometry. The z p hot technique is particularly useful to identify objects 
with similar SEDs, and thus to check the compatibility of multiple image con- 
figurations. A recent example has been published by Gavazzi et al. (2003), who 
derived the redshifts of the two main radial and tangential arc systems in the 
cluster MS2137-23 {z duster = 0.313) from UBVRIJK photometry, in order to 
improve the modeling of the central dark matter distribution. In this case, the 
redshift of the two arc systems is z p hot = 1.6 ± 0.1, in excellent agreement with 
the subsequent spectroscopic determination. z p h ot s can be used advantageously 
when photometric data span a large wavelength range, and particularly if near- 
IR data are available, because this allows to obtain accurate redshifts in the 
sensitive region of 0.8 < z < 2. 

Photometric redshifts are particularly useful to derive the mass from weak 
shear analysis. We can use z p hots to discriminate between cluster and foreground 
galaxies to scale the lensing distance modulus in order to compute the mass 
from the gravitational convergence. The average convergence k = S/S cr , which 
corresponds to the ratio between the surface mass density and the critical value 
for lensing, may be obtained as a function of the radial distance 9 using different 
methods (see Mellier 1999 for a review). The mass within an aperture is given 



where are the angular size distances between the cluster (I), the observer (o) 
and the source (s), and k(< 9) is the averaged convergence within the radius 9. 
Using the redshift distribution N(z p hot ) computed through a suitable filter set, 



of the mass. This method has been recently applied by Athreya et al. (2002) 
to determine the mass of the lensing cluster MS1008-1224. King et al. (2002) 
have used a similar procedure in A1689 to select background galaxies for their 
first detection of weak gravitational shear at infrared wavelengths. Because of 
the small projected surface across the redshift space, the effective surface which 
is "seen" through a cluster lens is relatively small, thus producing a strong 
variance from field to field. Obtaining the N(z p hot) distribution for each cluster 
helps improving the mass determination. 

Another different strategy for cluster mass reconstruction using lens magni- 
fication and z p hot techniques was proposed by Dye et al. (2001), with an applica- 
tion to Abell 1689. These authors used for the first time the lens magnification 
inferred from the luminosity function of background sources when scaling the 
convergence to real lens mass. 

Photometric redshifts are also particularly useful to improve the detection 
of clusters in wide-field surveys, as well as to properly identify the visible coun- 
terpart of complex and multiple lenses. Examples of composite lenses recently 
identified by z p hot are MS1008-1224 (Athreya et al. 2002), where a secondary 



by 



M(< 9) = k(< 9) S cr • vr (9D ol f = k{<9) 9 2 — 





can be computed, thus leading to a fair estimate 
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Figure 6. Comparison between spectroscopic and photometric red- 
shifts in the HDF fields (from Bolzonella et al. 2000). Photometric 
redshifts and 3er error bars were computed with hyperz. Most of the 
catastrophic identifications correspond to degenerate solutions. 
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lens exists at z ~ 0.9 in addition to the main cluster, and the Cloverleaf field 
(Kneib et al. 1998). Thompson et al. (2001) have used a deep optical and near- 
infrared survey to properly identify a massive cluster of galaxies at z= 1.263 with 
the X-ray source RX J105343+5735. Photometric redshifts techniques have been 
recently applied to the ESO Distant Cluster Survey (EDisCS, White et al. 2003, 
in preparation), which contains two different samples of clusters: a \ow-z sample 
at z ~ 0.5, for which BVIK photometry is available, and a high-z sample at 
z ~ 0.8, with VRIJK photometry. The first comparison between spectroscopic 
and photometric redshifts (on ~ 500 spectra) confirms the excellent behaviour 
of 

Zphot i the typical dispersion beeing u z ~ 0.06 and o z ~ 0.08 for the low and 
the high-z samples respectively. These results illustrate the typical accuracy 
that we can expect with well suited photometric data, at least up to a redshift 
of ~ 1. 

4. Looking for the first sources with Gravitational Telescopes 
4.1. Recent developments in PopIII studies 

One of the main possible applications of GTs is the research and study of the 
very first stars and galaxies forming from pristine matter in the early Universe, 
the so-called Population III objects (cf. review of Loeb & Barkana 2001; Weiss 
et al. 2000; Umemura & Susa 2001). These sources constitute the first building 
blocks of galaxies. Their detection still remains one of the major challenges 
of present day observational cosmology. The recent WMAP results seem to 
place these PopIII objects at redshifts up to z ~ 10-30 (Kogut et al. 2003), a 
redshift domain where the most relevant signatures are expected in the near-IR 
window. During the last years, modeling efforts have been motivated by future 
space facilities such as JWST, which should be able to observe these objects at 
redshifts z up to 30. Nevertheless, the detection and first studies on the physical 
properties of Population III objects could likely be started earlier using ground- 
based 10m class telescopes, and well suited near-IR facilities, such as Isaac/VLT, 
the multi-object spectrograph EMIR at GTC (~ 2006), or the future KMOS for 
the second generation of VLT instruments (~ 2008), at least up to z~ 18. 

According to modeling results, several observational signatures of Pop III 
stars and galaxies (i.e. ensembles/clusters of Pop III stars) are expected. In 
summary: 

• Strong UV emission and characteristic recombination lines of hydrogen 
and He II, especially Lyman a and HeII1640 (Tumlinson & Shull 2000, 
Bromm et al. 2001, Schaerer 2002, 2003). 

• Mid-IR molecular hydrogen lines at 2.12 /im and longer wavelengths formed 
in cooling shells (Ciardi & Ferrara 2001). 

• Individual supernovae whose visibility in the rest-frame optical and near- 
IR could be enhanced due to time dilatation (Miralda-Escude & Rees 1997, 
Heger et al. 2001). 

• High energy neutrinos from Pop III gamma-ray bursts possibly associated 
with fast X-ray transients (Schneider et al. 2002). 
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Rest-frame UV stellar and nebular continuous and recombination line emis- 
sion represent the largest fraction of the energy emitted by Population III ob- 
jects, which are generally thought to be predominantly massive or very massive 
stars, up to ~ IOOOA/q (e.g. Abel et al. 1998; Bromm et al. 2002; Nakamura & 
Umemura 2001; Nakamura & Umemura 2002). The predicted rest-frame UV to 
optical spectra of Pop III galaxies including the strongest emission lines (Lyman 
a, HeIIA1640 and others) have recently been computed by Schaerer (2002, 2003). 
The existence of a significant population of PopIII, or extremely metal-poor ob- 
jects, at redshifts z 10 is supported by several indications, including obser- 
vations of Ly-a emitters with unusually large equivalent widths in the LALA 
survey (Malhotra &: Rhoads 2002) and the recent simulations of Scannapieco et 
al. (2003) accounting for the non-homogeneous chemical evolution and pollution 
of the Universe. 

4.2. PopIII observational signatures 

Simulations have been carried out using the recent models by Schaerer (2002) to 
demonstrate the feasability of this scientific case in view of the future near-IR 
facilities mentioned above. The aim is to determine a well suited identification 
criteria for PopIII galaxies, and to study the feasability of spectroscopic obser- 
vations in the near-IR domain. Some preliminary results have been shown in 
Schaerer & Pello (2000) and Pello & Schaerer (2002). 

The SED of PopIII sources is dominated by the nebular continuous emission 
at A > 1400 A, and several strong emission lines are present, such as Lyman 
a, HeIIA1640, HeIIA3203, FMIA4686 and others. We have considered different 
IMFs and ages for the stellar population, as well as two different star formation 
regimes (starburst and continuous star formation), following the prescriptions 
given by Schaerer (2002). Lyman series troughs (Haiman & Loeb 1999), and 
Lyman forest following the prescription of Madau (1995) are included. The 
reionization redshift is assumed to be ~6. The virial radius is of the order of a 
few kpc for these sources, and thus we consider that sources are unresolved on 
a 0.3" scale, with spherical symetry. Different series of simulations have been 
performed to account for different hypothesis for the extended Ly a halo (cf. 
Loeb & Rybicki 1999). 

Figure Q displays the J-H versus H-Ks color-color diagram (Vega system) 
for Pop III objects obtained from these simulations, showing the expected po- 
sition of objects at different redshifts over the interval of z ~ 5 to 12, and the 
location for the most promising Pop III candidates. Various models for Pop III 
objects are presented in Figure [Tj In particular, models are computed assuming 
different fractions of Lyman a emission entering the integration aperture: 100% 
(thick solid lines), 50% (thin solid lines), and 0% (thin dotted lines). The ex- 
pected colors obtained within an integration aperture corresponding to the core 
(non-resolved) source, whith Lyman a emission coming from an extended halo 
(Loeb & Rybicki 1999), are displayed by thick dotted lines. The colors expected 
for a pure stellar population are given for comparison as dashed thick lines. We 
used a Salpeter IMF to derive the synthetic spectra in all cases. The location 
of stars and normal galaxies at < z < 8 is also given for comparison, as well 
as starbursts templates (SB1 and SB2, from Kinney et al. 1993, and the low 
metallicity galaxy SBS0335-052). The strong UV continuum and the presence 
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Figure 7. Top: J-H versus H-K' color-color diagram (Vega system) 
for Pop III candidates. Various models for Pop III objects are presented 
over the interval z ~ 5 to 12, with different fractions of Lyman a 
emission entering the integration aperture: 100% (thick solid lines), 
50% (thin solid lines), and 0% (thin dotted lines). Thick dotted lines 
correspond to colors obtained within an integration aperture limited to 
the non-resolved core of the source, and Lyman a emission coming from 
an extended halo (Loeb & Rybicki 1999). Positions of stars and normal 
galaxies up to z < 8 are also shown for comparison, including starbursts 
templates (SB1 and SB2, from Kinney et al. 1993 -blue dotted lines-, 
and the low metallicity galaxy SBS0335-052 -green dash-dotted lines-). 
The shift direction induced by reddening is also indicated. See text for 
more details. Bottom: Same color-color diagram as in Top Figure, 
showing the location of the iso-z lines: z < 8 (thick dashed lines), z=8 
(thin solid line), z=9 (thick dashed line), z=10 (thick solid line), z=10.5 
(thick dotted-dashed line) and z=ll (thin dashed line). 
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Figure 8. Expected number counts of PopIII objects with K'(Vega) 
< 24, in the field of view of Isaac, per 0.1 redshift bin interval, using 
different assumptions for the IMF (from Richard et al., in preparation). 
The typical results expected in the strong lensing regime (solid lines), 
towards the lensing cluster AC 114, are compared to the corresponding 
ones in blank fields surveys (dashed lines), for the same photometric 
depth. The gain in the z ~ 8 — 10 domain when taking advantage of 
the lensing amplification is typically a factor of 10 (see text for details). 
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of prominent emission lines allow to discriminate between PopIII objects with 
z > 8 and normal galaxies or stars. Even "normal" galaxies at z > 8 will exhibit 
redder H-K colors (which trace the restframe UV slope) as compared to PopIII. 
Because of the limiting magnitudes in the different filters, most PopIII candi- 
dates are expected to fall within the "PopIII candidates" region. This diagram 
illustrates that PopIII objects with strong ionising continua can well be sep- 
arated from "normal" objects provided that ultra-deep JHK photometry with 
sufficient accurary is available. 

We have also computed the expected S/N ratios for the main emission 
lines, as seen through a 10m telescope, using a medium resolution spectrograph 
(R~ 3000—5000). According to our results, Lyman a can easily be detected with 
a good S/N over the redshift intervals z ~ 8 to 18, with some gaps, depending on 
the spectral resolution (OH subtraction) and atmospheric transmission. A joint 
detection with Hell A1640, the strongest Hell line, is possible for z ~ 5.5-7.5 
(Lyman a in optical domain), z ~ 8-14 with both lines in the near-IR, again 
with some gaps. The typical line fluxes for the Hell A1640 line range between 
10~ 17 and a few 10 -18 erg/s/cm2, for a fiducial stellar mass halo of 10 7 Mq, 
corresponding to a collapsing dark matter halo of 2 x 10 8 Mq. The detection of 
both Hell A1640 and Ly a allows one e.g. to obtain a measure of the hardness 
of the ionising flux which constrains the upper end of the IMF and the age of 
Pop III systems. 

Thus, according to our results, efficient photometric and spectroscopic ob- 
servations of these features should be possible in the near-IR domain, thanks to 
the future imaging and spectroscopic facilities, for a large number of sources, 
thus allowing to derive statistically significant conclusions about their formation 
epoch and physical properties. Broad-band photometry is quite insensitive to 
changes in the IMF, and broad-band colors do not allow to constrain physical 
properties such as the IMF (i.e. the mass range of Pop III stars), ages, etc., 
but could be useful to identify the sources on ultra-deep photometric surveys. 
Spectroscopy is needed to study the physics of these objects. 

4.3. A search for PopIII candidates behind lensing clusters 

PopIII objects could be detected from deep near-IR photometry based on a 
measurement of two colors with accuracies of the order of 0.3 mag. Ideal 
fields for the first studies are lensing clusters of galaxies with areas of strong 
gravitational amplification. A devoted and prototype observational program is 
presently going on, aimed at the detection of Population III objects at high z 
using GTs. 

During the last year we were granted ISAAC /VLT time to obtain ultra-deep 
JHK' images of a gravitational lensing clusters in order to search for PopIII can- 
didates using our above mentioned color-color criteria. Two lensing clusters 
were observed up to now: AC114 (z = 0.31), images obtained late August 2002, 
and Abell 1835 (z = 0.253), observed in March 2003. At present, only data 
on AC 114 have been analyzed. Images were obtained under excellent seeing 
(0.4-0.6") and photometric conditions. The limiting magnitudes attained (3cr 
/ 4pixels) are J= 25.7, H= 24.7 and K'= 24.1 in the Vega system. From the 
first investigations on AC 114, several (~ 10) candidate PopIII galaxies are found 
in the selection window, with magnitudes K' ~ 23-24 (Vega system) and esti- 



Distant galaxies seen through Gravitational Telescopes 21 

mated redshifts of z ~ 8-11. All these objects are R-dropouts, located close 
to the theoretical position with respect to the critical lines. A few of them are 
brighter than H ~ 24 (including 1-2 mag amplification). They constitute the 
first spectroscopic targets for spectroscopic follow-up with ISAAC/VLT in the 
subsequent months. 

The expected number of Pop III objects and primordial QSOs in blank 
fields has been derived by several authors. E.g. in a comprehensive study Ciardi 
et al. (2000) show that at z > 8 naked stellar clusters, i.e. objects which have 
completely blown out their ISM, and thus avoid local chemical enrichment, dom- 
inate the population of luminous objects. Although pilot studies have recently 
started to explore the possible formation of dust in Pop III objects (Todini & 
Ferrara 2001), the effect is generally neglected. Based on such assumptions, Oh 
et al. (2001) have calculated the predicted number of Pop III objects detectable 
in Hell lines with JWST for a one day integration time. Their estimate yields 
between 60 and 4500 sources in a 10'xlO' field of view, depending on the model 
parameters. These estimates are still highly uncertain and model dependent. 

In Figure |S1 we compare the expected number counts of PopIII objects with 
K'< 24 in our GT ongoing survey, within the field of view of Isaac, per 0.1 
redshift bin interval, using different assumptions for the IMF. A simple Press- 
Schechter formalism for the abundance of halos and standard ACDM cosmology 
were considered in these simulations, as well as a conservative fixed fraction 
of the baryonic mass in halos converted into stars, assuming that all stars are 
PopIII, without chemical evolution (Richard et al., in preparation). Thus, num- 
ber counts shown in Figure E] have to be considered as a first order estimate, 
the important point beeing the comparison between lensing and blank fields. 
The typical results expected in the strong lensing regime, towards the lensing 
cluster AC114, are compared to the corresponding ones in blank fields surveys, 
for the same photometric depth. The gain in the z ~ 8 — 10 domain when taking 
advantage of the lensing amplification is typically a factor of 10. The expected 
density of primordial quasars could be similar to that of PopIII galaxies (Oh 
et al. 2001). Thus, multiplexing is needed in blank and lensing fields to allow 
highly efficient observations of relevant samples of Pop III objects. Also, because 
the observational signatures of primordial quasars are expected to be similar to 
those of genuine Pop III stars, a relatively high resolution is needed to obtain 
line profiles. 



5. Conclusions 

Gravitational Telescopes still remain a unique tool in deep universe studies. 
The recent developments in the sub-millimetric and mid-infrared domains, the 
search for extremely faint EROs and primordial galaxies, and the detailed study 
of the spectrophotometric and morphological properties of galaxies beyond the 
spectroscopic limits argue in favour of this technique. 

Recent results obtained by different authors on the physical properties of 
faint lensed galaxies (such as AC114-S2, AC114-A2, or the z ~ 3 sources in 
1E0657-558) suggest that high-z objects of different luminosities were subjected 
to different histories of star formation. But these results are based on a dramati- 
cally small number of objects. Enlarged samples of galaxies in both redshift and 
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luminosity, specially in the redshift range 1.5 < z < 6, are urgently needed to 
set strong constrains on galaxy evolution scenarios. The systematic use of GTs 
to access the most distant and faintest population of galaxies can help achieving 
this goal. 

An important issue for future spectroscopic studies of very high-z primordial 
galaxies is the strategy for target selection. Photometric redshifts and various 
color-selection techniques have demonstrated their ability to produce samples 
of high-z objects. Even young Pop III sources or extremely metal-deficient 
starbursts are expected to show distinct characteristics in their near-IR col- 
ors compared to "normal" galaxies at any redshift. GTs with areas of strong 
gravitational amplification are the ideal fields for the first prospective studies. 
Thanks to its multiplexing, spectral resolution and wide field-of-view capabil- 
ities, the forthcoming next generation of ground-based near-IR spectrographs, 
such as EMIR at Grantecan (http://www.iac.es/gtc/index.html) or the fu- 
ture KMOS for the VLT, are the suited instruments to complement JWST 
performances and to start exploring the formation epoch of the first stars in 
galaxies. 
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